Phase Separation and the Dual Nature of the Electronic Structure in Cuprates 
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The dual nature of the electronic structure of stripes in La2- x Sr x Cu04 was characterized by many 
experiments. We present here an attempt to characterize this dual behavior based on the Cahn- 
Hilliard theory of a phase separation transition which is assumed to occur at the high pseudogap 
temperature. The resulting inhomogeneous low doping system is formed of hole-rich (metallic like) 
regions embedded in a hole-poor (insulator like). This inhomogeneous configuration is analyzed by 
a new method within the Bogoliubov-deGennes superconducting theory. This approach describes 
well the electronic nodal-antinodal dichotomy and parts of the phase diagram. 

PACS numbers: 74.72.-h, 74.80.-g, 74.20.-z, 02.70.Bf 



The high critical temperature superconductors 
(HTSC) represents today one of the greatest challenge 
of condensed matter physics. It is likely that the main 
difficulty to understand the physical properties of HTSC 
is due to the fact that some families seem to have a high 
inhomogeneous electronic structure while others families 
appear to be more homogeneous or, at least, without 
some gross inhomogeneity. 

An important technique to study the HTSC electronic 
structure is provided by angle resolved photoemission 
(ARPES) experiments. With the improvement of the 
energy and momentum resolution in recent works 0, U 
Eli , it was possible to distinguish a two component 
electronic structure in fc-space of the Lai- x Sr x CuO± 
(LSCO) family: a metallic quasi particle peak crossing 
the Fermi level along the zone diagonal following the 
(0, 0)-(tt, 7t) nodal direction which increases with the dop- 
ing levelUI^- On the other hand, the spectral weight 
at the straight segments in the (w, 0) and (0,7r) antin- 
odal regions, is compatible with a quasi one dimensional 
structure or with static stripes^. Due to the d-wave 
symmetry of the superconducting order parameter, the 
zero temperature superconducting gap Ao vanishes at 
the nodal and is measured at the antinodal directions 
by the leading edge shift on the Fermi surface by the 
ARPES spectra. Moreover, the values of A decreases 
with doping but the quasiparticle spectral weight near 
the nodal directions, at the Fermi level, increases, show- 
ing the distinct behavior of these two aspects of the elec- 
tronic structured H 0, || . 

A complementary technique to ARPES is provided by 
scanning tunneling microscopy (STM) since it probes the 
differential conductance or Ao directly on the surface of 
the compound. Recent STM data have revealed a patch- 
work of (nanoscale) local spatial variations in the density 
of states which is used to measure the local supercon- 
ducting gap0,@- With this technique, it was also pos- 
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sible to distinguish two distinct behavior: well defined 
coherent and ill-defined incoherent peaks depending on 
the spectra location on a Bi2Sr2CaCu20s+s (Bi2212) 
surface |9l lid. Hl|. Also, tunneling experiments using su- 
perconductor insulator superconductor (SIS) with differ- 
ent insulator layers have shown distinct sets of energy 
scales and have also led to the idea that the richness of 
the phase diagram as function of dopin g is due to the 
charge inhomogeneity in the Cu-0 planes 12]. 

These and others unusual features of cuprates led to 
theoretical proposals that phase separation is essential 
to understand their physics [13, 13 ■ I n f ac ti phase segre- 
gation has been observed on the Lo^CilOa+s system by 
x-ray and transport measurements [l5|, [lfj. They have 
measured a spinodal phase segregation into an oxygen- 
rich (and hole-rich) metallic phase and an oxygen-poor 
antiferromagnetic phase above T=220K. Below this tem- 
perature the mobility of the interstitial oxygen becomes 
too low for a further segregation. La^CuO^+g is the only 
system where ion diffusion has been firmly established, 
although there are evidence of ion diffusion at room tem- 
perature in micro crystals of the Bi2212 superconductors 
at a very slow rateflTj. 

In this paper we take the large pseudogap tempera- 
ture of the HTSC phase diagrampJI [ljj as the phase 
separation temperature T ps (p). To obtain quantitative 
results on the charge separation as a function of the tem- 
perature, we apply the Cahn-Hilliard (CH) theory 20] to 
this transition. It yields the phase separation patterns 
found in many HTSC and it provides also an interpre- 
tation to the energy E g , common to all cuprates |is|. as 
the potential barrier between the two equilibrium hole- 
rich and hole-poor phases. Since T ps (p) increases as the 
average doping level p decreases, the differences between 
high and low local densities are enhanced in low doping 
compounds. With the charge structure that comes out 
of these calculations in a N x N cluster, we use the BdG 
method, with the local chemical potential associated with 
this variable charge structure, to calculate the supercon- 
ducting properties. This is a novel approach to the charge 
inhomogeneities in the BdG context. Several interesting 
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features comes out of this procedure which allows us to 
study the insulator, metallic and superconducting phases 
by quantitative calculations that compares well with the 
experimental results. 

The CH non-linear differential equation which de- 
scribes the time evolution process of a phase separation 
process, at a temperature T below the phase separa- 
tion transition at T ps (p), can be written as a time (t) 



derivative [191 12 



^ = -AfV 2 (e 2 W 
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where u is the local order parameter. In a finite size 
scheme pH|, it is associated with the local charge p(i), i.e., 
the local variation from the average number of holes per 
copper atom p, defined at a site i as u(i) = p(i) — p and 
u(i) « above and near the T ps . e = 0.01 and B — 1 are 
fixed parameters, A depends on the temperature T and 
the ratio ±A/ B yields the two local equilibrium densities 
p± , a low and a high density. M is the mobility of the 
particles and it dictates the phase separation time scale. 
Thus, the degree of phase separation depends on how 
fast the system is quenched, what is an explanation why 
similar compounds may exhibit different degree of inho- 
mogeneity. As the temperature goes down below T ps , the 
two equilibrium order parameters (or densities) increases 
their differences and the energy barrier between them, 
Eg, also increases E g 




A (T)/B which is propor- 
tional to (T ps — T) 2 and can be taken as the energy as- 
sociated with the upper pseudogap temperature [l9|. In 
Fig. we display the mapping of the order parameter for 
a 100 x 100 system with p = 1/8. The two equilibrium lo- 
cal densities for this system is p- = zero and p + = 0.24. 
At the beginning, the phase separation process is quite 
symmetric, and could be the reason for the symmet- 
ric phases as stripes and checkerboard |22^ order. As it 
evolves in time, the systems tends towards a complete 
phase separation, with larger and less ordered stripes, as 
shown in Fig.QJ. 

Our main point here is that this phase separation pro- 
cess is the cause of the charge inhomogeneity found in 
many HTSC materials. Taking the large pseudogap, as 
the phase separation temperature T ps (p), it is possible 
to infer that a HTSC compound, depending on its dop- 
ing level, evolves in different patterns as, for instance, 
in the way depicted in Fig.QJ. Therefore underdoped 
compounds which, in general, have a very high T ps , may 
phase separate into a complete bimodal charge distribu- 
tion with larger charge stripes or phase domains (Fig. 
or d)), while compounds with j)« 0.19 which have very 
low values of T ps , segregates in small and symmetric re- tion) operators at site Xj, spin a{] J,}, and = c\ a Ci a . 
gions as shown in Fig.Qi or b). Above T ps or for com- t%j is the hopping between site i and j. In an at- 



FIG. 1: (Color online) The mapping of the order parameter 
during the process of phase separation in a compound with 
average doping p = 0.12, starting with variations around u = 
0. Panel a) is for t=400 time steps, b) t=800, c) t=1000 and 
d) t=4000. 



perature pseudogap peaks (ZTPG) for compounds ap- 
proaching 0.19. In the CH scenario, the origin of the 
ZTPG is the value of the ener gy b arrier E g between the 
hole-poor and hole- rich phases[l!j- Since the charge do- 
mains and Eg decrease with p, the ZTPG must behave 
in a similar fashion: decreasing in number and intensity 
and vanishing near p = 0.19, exactly as observed in re- 
cent STM resultsHflllJ. 

Now, that we have derived how the charge inhomogene- 
ity sets in a HTSC material below the high pseudogap 
temperature T ps , we can study how the superconduc- 
tivity develops in a charge inhomogeneous system, like 
in a stripe phase, as function of the temperature. With 
this aim, we feed into a local superconducting calculation 
with the BdG mean- field theoryHH, the charge or do- 
mains mappings derived from the CH solutions. As usual, 
the method starts with the extended Hubbard Hamilto- 



h = - % c L c jv + J2^ 7 
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where c\ a (c icJ ) is the usual fermionic creation (annihila- 



pounds with p > 0.19 the systems are described by small 
doping variations around the mean value or a Gaussian 
distribution. This change in behavior near p sa 0.19 has 
been detected in many experiments^^]. More recently, 
STM data|T3,II3 detected the vanishing of the zero tem- 



tempt to model real systems, we have used hopping val- 
ues up to 5 th neighbors derived from the ARPES data 
of YBCOp5l|. In their notation, the hopping parame- 
ters here are: t = ti=0.15eV, t 2 /*i=-0.70, i 3 /*i=0.25, 
i 4 /ti=0.08, t 5 /ti=-0.08. U = Lit is the on-site and 
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V = -l.lt is the nearest neighbor phenomenological 
interactions, fa is the local variable chemical poten- 
tial which reproduces the charge inhomogeneity solutions 
shown in Fig.©, according the CH results. The BdG 
equations are: 



(3) 



These equations are solved self-consistently in clusters 
from 14 x 14 to 24 x 24 sites for the positive eigenvalues 
E n or quasiparticle excitations and the eigenvectors or 
charge amplitudes UjApCj ) and w n (xi), together with the 
pairing amplitudes |23j 



t=0.lSeV t 2 =-0.70t t j=0.2St t 4 =0. 




Aj/(xi) = -U^2u n (-Xi)v*(-Xi) tanh 



EL 
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V \ 

A (5 (x i ) = -— 2^K(xi)<(xi +5) 



+v* (x i )u„(x i + J)] tanh 
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and the hole density is given by 



JK(x*)| 2 /« + Mx 4 )| 2 (l -fn)}, (6) 



where /„ is the Fermi function. 

The major difference from previous calculations is in 
the way that the disorder is taken into account. Ghosal 
et al|24j used an impurity potential defined by a random 
variable between the limits [— V, V]. Here, as mentioned 
above, we feed into the initial conditions the non-constant 
local charge, as those shown in Fig.©, derived from the 
CH solutions. With this fixed charge mapping in a N x 
N cluster, we let the chemical potential fa change, self- 
consistently, until it yields these initial charge domains. 

Fig-© shows the superconducting d-wave gap A(i,T) 
at each site i of two 14x21 clusters and their temperature 
evolution. Panel a) represents an underdoped compound 
of p = 0.05 with total phase separation into a bimodal 
charge distribution made of charge variable stripes of 14 
sites each, and in a less symmetric geometry, as taken 
from Fig.©i). The values of the local doping level p(i) 
are shown on the top of each panel, p = 0.0 is the light 
color and p(i) = 0.24 is the dark color phase depicted in 
the Fig.©i). Panel b) represents a compound of average 
doping of p = 0.12, made of values of p(i) = 0.0, 0.12 and 
0.24 similar to Fig.©). The gaps are calculated by the 
BdG equations and the chemical potential p,(i) evolves 
self-consistently at each site i, until it yields the initial 
conditions on the local densities, which are held fixed. 
We have verified that changing the stripes configuration 
in the clusters yields the same type of gap structure and 
temperature dependence. 
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FIG. 2: (Color online) The temperature evolution of the 
local gap A(i,T) (in units of t = 0.15eV) with stripe charge 
disorder (the values at each site i are given on the top of 
each panel) and average doping level p = 0.05 and 0.12. The 
horizontal axis displays 21 stripes with 14 unit cells each. 



The results for the p = 0.05 sample with p_ = and 
p + = 0.24 is shown in the top panel of Fig.©. As one 
can see, the largest values of the zero temperature gap 
Ao(z) is in the regions where p(i) = 0.24. This is a di- 
rect interpretation of the ID metallic behavior and the 
corresponding high values of the leading edge shift on 
the Fermi surfaceQ along the (jr, 0) and (0, tt) antinodal 
straight segments. As p increases, the size of hole-rich 
stripes increases, for instance, for p = 0.05 they are typi- 
cally made of 3 lines and for p = 0.12 they are made of 6 
lines. Thus, upon doping, the hole-rich stripes increases 
in number and size, changing the properties of the system 
from one dimensional to two dimensional character and 
enhances the overall metallic electronic behavior. Such 
changing with the average doping level was also verified 
by the ARPES data on many lightly doping samples 2] by 
the measurement of spectral weight at the antinodal and, 
the increase of the spectral weight along the (±tt, ±7t) 
nodal regions with p. 

We discuss now how T c (p) and T*(p) (the lower pseu- 
dogap) are estimated in our calculations. We see from 
Fig.© that, at zero temperature, there is a large vari- 
ation of A(i,T) throughout the samples. Hole-poor re- 
gions have A(i,T) « and hole-rich have finite values. 
As the temperature increases slightly, mostly local gaps 
at the hole-poor regions vanish, leaving only finite su- 
perconducting gaps at the hole rich regions. Thus the 
p = 0.05 sample, as concerns it resistivity, is an insu- 
lator at all finite temperatures, although it has a few 
metallic stripes which were detected by ARPES [2|- For 
the p = 0.12 compound, even the gaps at the hole-poor 
stripes remain up to T = 40K. Thus, it is clear that it is 
the superconducting percolating temperature which, in 
agreement with previous work|26|. is the superconduct- 
ing critical temperature T c . In other words, above T c , the 
local superconducting gaps vanish at hole-poor regions 
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and it prevents the superconducting regions to percolate, 
there are still some non- vanishing A(i,T) at the hole- rich 
or metallic regions but, since they occupy marginally less 
than 50% of the system size, it cannot hold a supercon- 
ducting current. Thus, above T c , the p = 0.12 sample 
has a metallic behavior. As it is also shown in Fig.J2J), 
these local gaps decreases continuously as the tempera- 
ture increases further and, for p = 0.05 and p = 0.12, 
they totally vanish at T = 110K and T = 90K respec- 
tively. This is the interpretation to the onset of the local 
superconducting or lower pseudogap temperature T*{p). 
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FIG. 3: (Color online) The maximum local zero temperature 
gap Ao(i) (squares) compared with the experimental leading 
edge shifty (circles) as function of p. It is also shown the 
phase separation line or upper pseudogap energy (E g ). In the 
inset A(T) for the stripe and homogeneous system. All the 
gaps and energies are in units of t — 0.15eV . 

The maximum zero temperature Ao(i) as function of 
p is shown in Fig.©- Aq(«) scales reasonable well with 



the ARPES measurements of the leading edge shift Q 
demonstrating that the maximum gap in a sample in- 
creases with the disorder. As already shown in Fig.©, a 
system with density disorder in the form of stripes, en- 
hances the local gap. Thus, we find that the maximum 
gap for a such inhomogeneous sample is larger than that 
of a homogeneous similar (same p), that is, the disorder 
enhances the local gaps at the hole-rich regions and at 
their borders. However, with the same coupling, the av- 
erage gap in a disordered system is less than the gap of 
a homogeneous one, and the average values of T c is also 
enhanced. These results are shown, for s-wave calcula- 
tions, in the inset of Fig.Q and are in agreement with 
some recent calculations 27]. 

In summary we have used the CH phase separation 
approach to model the inhomogeneity found in LSCO. 
The derived charge stripe-like disorder describes well the 
dual nature of the underdoped electronic structure. The 
energy barrier E g between the two equilibrium densities 
furnishes an interpretation to the large pseudogap energy 
scale and the ZTPG measured on Bi2212 by STM. The 
spinodal separation provides also a mechanism of how 
the electronic structure of LSCO, as seen by ARPES, 
evolves with doping; the increase of the leading edge shift 
at the straight (ID) segments near the antinodal region 
and the increase of the spectral weight near the nodal 
region, with the concomitant developing of the metallic 
behavior. The disorder favor the clustering of metallic re- 
gions embedded in an insulator matrix at the low doping 
compounds. The site dependent local superconducting 
gaps at low and high doping regions provide also an in- 
terpretation to the onset of superconductivity or lower 
pseudogap T*(p), in agreement with the phase diagram 
and the nonconventional properties of HTSC. 
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CNPq and CNPq-Faperj Pronex E-26/171. 168/2003. 
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